Variational Approach to the Spin-boson Model With a Sub-Ohmic Bath 
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The influence of dissipation on quantum tunneling in the spin-boson model with a sub-Ohmic bath 
is studied by a variational calculation. By examining the evolution of solutions of the variational 
equation with the coupling strength near the phase boundary, we are able to present a scenario 
of discontinuous transition in sub-Ohmic dissipation case in accord with Ginzburg-Landau theory. 
Based on the constructed picture, it is shown that the critical point found in the general way is not 
thermodynamically the critical point, but the point where the second energy minimum begins to 
develop. The true cross-over point is calculated and the obtained phase diagram is in agreement 
with the result of numerical renormalization group calculation. 
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I. INTRODUCTION 



The spin-boson model is an important toy model for 
investigating the influence of dissipation on quantum 
tunneling and has a wide range of applications. Q 
Over decades the model has been studied by various 
methods, like path integral, renormalization group 
calculation, [4, la variational calculation, 
[m IT^ and the numerical renormalization group(NRG) 
calculation. [l^j etc.. One important issue is to study 
the cross-over from the delocalized to localized phases 
as the dissipation increases. Most of the studies are con- 
centrated on the Ohmic dissipation case which is con- 
sidered as corresponding to real physical systems and 
the cross-over picture is well understood. On the other 
hand, the situation for the sub-Ohmic dissipation, which 
is of less physics interest but still important for a well 
understanding of the spin-boson model, has some con- 
fusions. Renormalization group calculation shown that 
quantum tunneling is totally suppressed by dissipation 
for any non-zero sub-Ohmic coupling at T=O,|3,0 while 
different conclusion was found by mapped the spin-boson 
model to an Ising model [l^ and using the well-known 
result for Ising model. The sub-Ohmic case was also 
studied by using infinitesimal unitary transformation and 
the cross-over was found to be discontinuous. Re- 
cently, the NRG calculation, which is considered as a 
powerful tool for investigation of the Kondo model and 
its generalizations, confirmed the delocalized to localized 
phase cross-over in sub-Ohmic dissipation case and the 
cross-over is identified as continuous. Variational cal- 
culation has been used to study the spin-boson model 
with a Ohmic bath and the result of cross-over bound- 
ary is in good agreement with the renormalization group 
calculation. 0, ^ The variational calculation for 
non-zero temperature j9j was generalized to sub-Ohmic 
case recently and the discontinuous cross-over behavior 
was found to exist at non-zero-temperature. J.2j Up to 
now, the description for this discontinuous cross-over is 
just limited to the discontinuous change of the tunnel- 
ing splitting at the cross-over point, while a scenario for 
such a discontinuous cross-over is still lacking. Accord- 



ing to Ginzburg-Landau theory It', IS^, the evolution of 
the free energy around the critical point for the first or- 
der (discontinuous) phase transition is rather complicated 
and merely a discontinuous change of order parameter at 
the cross-over point is certainly no enough for a com- 
plete description of this discontinuous transition. In this 
paper, we present further analysis on this discontinuous 
cross-over by examining the evolution of the solutions of 
the self-consistent equation derived from the variational 
calculation. It is found that the evolution of the solutions 
near the phase boundary is consistent with the general 
picture of the first order phase transition. Basing on the 
constructed picture, it is shown that the critical points 
determined in the general way are not thermodynami- 
cally critical points and the true critical point is calcu- 
lated. The arrangement of the paper is as follows. In 
the next section, the model and a brief explanation on 
variational calculation are presented. In section III we 
present analysis on the discontinuous phase transition 
by comparing the evolution of the solutions of the self- 
consistent equation for Ohmic and sub-Ohmic dissipation 
cases near the critical point. Conclusions and discussion 
are given in the last section. 



II. THE MODEL AND VARIATIONAL 
CALCULATION 



The Hamiltonian of the spin-boson model is given 
by (setting n= 1) 0,12 

H = |(T^ -t- ^(Tx + ^ h\hkuik + ^ Ck{h\ + bk), (1) 



where (Ji{i = x,y,z) is the Pauli matrix, bk{bj.) is the 
annihilation(creation) operator of the kth phonon mode 
with energy Uk and Ck is the coupling parameter. The 
main interest of the present paper will be the zero tem- 
perature so we set the bias e = in the following. It is 
known that the solution of this model is determined by 
the so-called the bath spectral function(density) defined 
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JH = ^^45(c.-c.fe). (2) 

k 

Generally J{uj) is characterized by a cut-off frequency lOc 
and has a power-law form, i.e., 

J{lu) = |awVw^-\ < w < cjc, (3) 

where a is a dimensionless coupling strength which char- 
acterizes the dissipation strength. Parameter s specifies 
the property of the bath, s = 1 is the case of Ohmic 
dissipation and < s < 1 the sub-Ohmic dissipation 
case. It should be noted that J(w) can take some dif- 
ferent forms,[l[l3 like = f a^Vw^-ie""/'^- and 
J2{^) = ^auj" /ujg~^ with ojc — > c», while we find that 
the solution is almost the same in sub-Ohmic case(see 
below). 

As one can see from the Hamiltonian given in Eq.(l), 
when A = the localized phase is stable since in this case 
we have [<Jz, H] = 0. This result implies that the coupling 
to phonon bath alone cannot lead to tunneling and thus 
this problem can be treated approximately in the way 
without coupling to the bath as given in the quantum 
mechanics textbook. 19] To ensure the tunneling is small 
which is a precondition of our treatment, the following 
calculation is restricted to the condition A/wc <C 1. We 
denote the eigen-states of spin-up(down)-plus-bath as | t 
)|(/)+)(| where \(f)±) represent the eigen-state of 

the phonon bath without tunneling. When the tunneling 
is taken into account, the eigen-state of the whole system 
can be approximately given bv|l9| 

|<i>±) = (|T)|0+)±U)|0-))/V2, (4) 

then the tunneling splitting in the presence of dissipation 
is 

A' = ($+11/1$+) - = A(0+|(/._), (5) 

a well known result that A' is determined by the over- 
lap integral of the phonon ground states. [3. 0| In 
the absence of tunneling(i.e., A = 0), Hamiltonian (1) 
can be diagonalized by a well-known displaced-oscillator- 
transformation and the phonon ground states are the so- 
called displaced-oscillator-statespoj 

|</),±)=exp{±^^(&.-5l)}|0), 

where |0) is the vacuum state of the phonon. The essence 
of the variational calculation is that, in the presence of 
tunneling, the phonon ground states are suggested to still 
have the same from, i.e., 

|0±) =exp{±^gfe(6fe-6l)}|0), (6) 

k 

but leaving the parameter to be determined from the 
condition that the ground state energy of the whole sys- 
tem is a minimum with respect to gk- Substituting the 



above equation to Eq. (4) , the ground state energy of the 
whole system is found to be 

E[9k] = - 2cfcgfe) - ^A exp{-2 ^ 5^1^ (7) 

k k 

which is a functional of g^ , then = leads to 



^'^ c., + Acxp{-2Efe5^}' 
the tunneling splitting, by Eq.(5), is given by 

A' = ifA, X = F[.gfe] = exp{-2^.g2}. (9) 

k 

Using Eq.(8) and the definition of the spectral func- 
tion, we find that K is determined by the following self- 
consistent (or variational) equation 

K = f{K), f{K) ^ exp{-a /' ^fj\,^,, }- 

(10) 

Such a kind of self-consistent equation has been derived 
in previous works 0, H, S Efl and it plays an im- 

portant role in dealing with the cross-over from the delo- 
calized to localized phases. It is easy to see that K = 
is the trivial solution of Eq.(lO) and this solution repre- 
sents the localized phase. When the coupling strength a 
is large enough, K = Q is the only solution of the self- 
consistent equation, while as a decreases to some value 
ttcj the self-consistent equation begins to have, in addi- 
tion to the trivial solution, non-zero solutions, then etc 
is identified as the critical point at where the cross-over 
from the localized(i^r = 0) to delocalized(-ftr > 0) phases 
happens. This is the general way to determine the phase 
boundary used in previous works. [sl I9I ITol fill IT^ 

In the case of Ohmic dissipation, i.e., s = 1, the self- 
consistent equation can be solved analytically and the 
phase boundary is found to be ac = 1 in the case of 
A/ LOc <C 1, in agreement with the renormalization group 
calculation. f9lfTl| In the case of sub-Ohmic dissipation, 
the self-consistent equation can be solved numerically 
and the phase boundary between the localized(iir = 0) 
and delocalized(A' > 0) phases determined in this way 
is shown in Fig.l. The result by using different spectral 
functions, i.e., Ji((jj) and J2(w), are also shown in the 
inset. Our result shows that, for A/wc ^ 1, the phase 
boundary is almost the same for all three spectral func- 
tions as s < 0.7, while Uc is a little bit lower for J2ii^) 
when s > 0.7. Also, it is found that the relation between 
the critical coupling ac and A /ujc has a simple power-law 
form ac cx (A/wc)^"" as found by NRG calculation. [T^ 
Notably, such a relation can be deduced from Eq.(16) 
in ref . |la| by using the spectral function given here.|2l| 
However, as we shall show in the next section, the ac 
determined in this way for sub-Ohmic case is not ther- 
modynamically the critical point, but just the the limit 
of metastability for superheating of the first order phase 
transition. 1^] 
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FIG. 1: Phase boundary determined by for various A/ojc- 
The inset shows the comparison with the result by using dif- 
ferent spectral function Ji{ui) and J2(a;)(see the text) in the 
case of A/ujc — 0.01. 
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FIG. 2: Evolution of the solutions of Eq.(lO) with the increase 
of coupling strength a for A/ujc — 0.1 in the case of Ohmic 
dissipation s — 1. 



III. THE DISCONTINUOUS CROSS-OVER IN 
SUB-OHMIC CASE 

Now we turn to present a scenario for such a dis- 
continuous cross-over from the delocalized to localized 
phases in sub-Ohmic case. The key point is to exam- 
ine the evolution of the solutions of the self-consistent 
equation with the coupling strength a near the phase 
boundary. For clarity, we first see what happens in the 
Ohmic case. Fig. 2 shows the evolution of the solutions of 
Eq.(lO) with the increase of a in Ohmic dissipation case. 
When a > ttc, we have the trivial solution only, while a 
non-zero solution(/iri ^ 0) appears for a < ac- As one 
can see from the figure, the non-zero solution Ki contin- 
uously tends to as a approaches ac- This is consistent 
with the picture of a continuous (second order) transition: 
[TtL IT^ above the critical point (a > Uc), there is only 
one stable phase(one energy minimum located at some 
gko satisfying F[gko] = in the present case), below the 
critical point, this stable phase becomes unstable(i?[gfco] 
becomes the maximum of the energy) and a second stable 
phase appears(£^[(7fei] is the new energy minimum, where 
-P'lSfei] = Ki > 0), the cross-over behavior is continuous. 

The situation for the sub-Ohmic dissipation case is 
qualitatively different. As shown in Fig. 3, when a < ac, 
there are tiiio non-zero solutions of Eq.(10)(i4r2 > Ki 0) 
in additional to the trivial solution. 22] As the coupling 
strength a increases, K2 decreases while Ki increases and 
tends to meet K2 as a approaches ac- At a = ac — 0, we 
have Ki — K2 = Kq ^ and at this point Kq is the point 
of tangency for the line y — x and curve y — f{x)- At 
a — ac-\-Q, the solution Kq disappears suddenly and only 
the trivial solution is found. The a-dependence of the 
non-zero solutions of Eq.(lO) for Ohmic and sub-Ohmic 
dissipation cases are shown in Fig. 4 where one can see 
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FIG. 3: Evolution of the solutions of Eq. (10) with the increase 
of coupling strength a for A/uc — 0.1 in the case of sub- 
Ohmic dissipation s = 0.3. 
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FIG. 4: The a-dependence of the non-zero solution of Eq.(lO) 
for s — l(left) and s — 0.3(right). As a/ac — > 1, the non-zero 
solution of s = 1 approaches continuously, while for s = 0.3, 
the non-zero solution jumps from Ko 7^ to 0. 
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FIG. 5: Evolution of the K-dependence of the ground state 
energy E with coupling strength a for sub-Ohmic case(s = 
0.3) when A/ujc = 0.1. 



FIG. 7: The phase boundary determined by etc and ai in the 
case of A/lOc = 0.1. 




F\<,k2] = K2 



FIG. 6: The evolution of the energy extrema with the coupling 
strength a for the sub-Ohmic case. Thermodynamically the 
cross-over point is qi, while Qc is just the point where the 
second minimum begins to develop. 



the qualitatively different cross-over behavior. The re- 
sult of sub-Ohmic dissipation case clearly shows that the 
cross-over is discontinuous since the non-zero solution of 
Eq.(lO) and thus the tunneling splitting A' changes dis- 
continuously at the point a = ac- Such a behavior was 
found before [H [p and took as the evidence for a dis- 
continuous transition since the tunneling splitting has a 
physics meaning of the order parameter. 

What we want to emphasize here is the two non-zero 
solutions when a < ac- Physically we need to know 
which solution is stable and the meaning of the second 
non-zero solution. As one can see from Fig.4(right), 
the non-zero solution Ki increases with a, one can in- 
tuitively conclude that Ki is unstable since physically 
the tunneling splitting should decrease with a. Bas- 



ing on the variational calculation, we cannot make fur- 
ther analysis on the stability of the solution, so we turn 
to energy analysis, i.e., to see which solution is energy 
preferable. Typical evolution of the E — K curve with 
a is shown in Fig. 5. The result shows that, we have 
E{Ki) > E{0) and E{Ki) > E{K2) when a < 
E{0) decreases while both E{Ki) and E[K2) increase 
relatively as a increases but E{Ki) is always the largest, 
finally E{Ki) = E{K2) = E{Kq) as a = ac and we 
have £;(0) < E{Kq). This implies that, both K = Q 
and K2 are energy preferable while Ki is unstable when 
a < ac- Such a result is consistent with the scenario of 
a first order transition. In the scenario of the first order 
transition,'!^, Ts'l below the critical point, there are two 
free energy minima and a maximum lies between, while 
above the critical point, only one global free energy min- 
imum survives. In the sub-Ohmic dissipation case, when 
a < ac, three solutions of Eq.(lO) represent the two en- 
ergy minima and one energy maximum, that is, E[gko] 
with F[gko] = and E[gk2] with F[gk2] = K2 are the 
two energy minima, while i?[5fci] with F[gki] — Ki is the 
energy maximum lies between as shown in Fig. 6. As a 
increases and approaches a^ E[gki] tends to meet E[gk2] 
and at the point a — ac, these two energy extrema merge 
into a point of inflection at -^[(7*;] = Kq, then only one 
energy minimum £^[17^0] survives when a > ac- Based 
on the picture for the discontinuous phase transition, 
it is now clear that ac is not the critical point for the 
cross-over to happen, but just the point where the sec- 
ond energy minimum begins to develop, ac can be con- 
sidered as the limit of metastability for superheating. [Tsj 
i.e., the limit of metastability for increasing the dissipa- 
tion strength in the present case. Thermodynamically 
the critical point, as shown in Fig. 6, should be ai where 
we have m 13 



E{K2) = E{Q), 



(11) 
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FIG. 8: The phase boundary determined by a\ for various 
A/oJc- The inset shows the comparison with the NRG calcu- 
lation in the case of A/ci^c = 0.001. 



from this, (X\ can be determined by Eqs.(7), (8) and (10). 
Comparison between the phase boundary determined by 
ttc and a\ is shown in Fig. 7. It is easy to see that a\ < ac 
while the difference between ac and ai decreases as s in- 
creases and tends to zero as s ^ 1 where the transition 
becomes continuous. We also find that the difference be- 
tween ac and ai decreases with A/lUc- The phase bound- 
ary deduced in this way is shown in Fig. 8 which is simi- 
lar to that shown in Fig.l but with all the critical points 
lower. It is found that the phase boundary determined 
by ai is in good agreement with that obtained by NRG 
calculation when A/ujc < 0.01. 



IV. CONCLUSIONS AND DISCUSSION 

In conclusion, we have study the cross-over behav- 
ior from localized to delocalized phases of a spin-boson 
model with a sub-Ohmic bath by variational method. By 
examining the evolution of the solutions of self-consistent 
equation (10) with the coupling strength, we are able to 
present the scenario of the discontinuous transition in 
sub-Ohmic dissipation case. Based on the constructed 
picture, it is shown that the ac, at where the self- 
consistent equation begins to have non-zero solutions, is 
not thermodynamically the critical point, but just the 



point where the second energy minimum begins to de- 
velop. The true critical point is determined according to 
Ginzburg-Landau theory for the first order phase tran- 
sition and the obtained phase boundary is in agreement 
with the NRG calculation. Our analysis shows that the 
cross-over behavior in spin-boson model is directly re- 
lated to the evolution of solutions of the self-consistent 
equation derived from the variational calculation. The 
evolution behavior of solutions for a continuous cross- 
over (in Ohmic dissipation case) is qualitatively different 
from that of a discontinuous cross-over(in sub-Ohmic dis- 
sipation case). The present work, on one hand, provides 
convincing evidence for a discontinuous cross-over in sub- 
Ohmic case and on the other hand, demonstrates the new 
way to deal with the cross-over behavior in spin-boson 
model by the variational method. 

According to the definition of stable and unstable fixed 
points for renormalization group, 23] geometrically one 
can see from Fig. 3 that, both K — and K2 are stable 
fixed points while Ki is unstable fixed point as a < ac 
in sub-Ohmic case. On the other hand, we only have one 
stable fixed point(i.e., Ki) and one unstable fixed point 
as a < ac in Ohmic case. This result is in agreement with 
the NRG calculation, where 3 fixed points (2 stable 
and 1 unstable) were found in sub-Ohmic case while the 
third unstable fixed point disappeared in Ohmic case. 
However, the cross-over behavior in sub-Ohmic case was 
identified as continuous in NRG calculation, this implies 
further analysis is needed for seeking a deeper relation. 
Although the work by Kehrein and Mielke is not based 
on the variational calculation,^^ the cross-over behavior 
was studied by a self-consistent equation and the discon- 
tinuous behavior was judged by the discontinuous change 
of the tunneling splitting at the critical point ac, where 
the self-consistent equation begins to have non-zero solu- 
tions. Some results, like the (A/utc) dependence of criti- 
cal coupling ac and the s-dependence of tunneling split- 
ting at the critical point also show quantitative agree- 
ment with our work determined from Eq.(lO) at a — ac- 
This may lead to a conclusion the the critical point de- 
termined in ref.|16j is just ac given in the present work, 
i.e., not thermodynamically the critical point. 
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